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Summary 


A theoretical method is presented for predicting the lateral- 
directional stability derivatives of wing-body combinations with or 
without the blowing jet effect. The fuselage effect is accounted 
for by the axial distribution of G.N. Ward's vortex multiplets. 
Comparison of the predicted results with experiment and other 
theoretical methods shows good agreement for configurations without 
the blowing jet. More applicable experimental data with blowing 
jets are needed to establish the accuracy of the theory. 



1. List of Symbols 


AR 

aspect ratio 

b 

span 

c 

chord length 

C DJ 

induced drag coefficient 

C L 

lift coefficient 

C * 

rolling moment coefficient 


pitching moment coefficient about 

C N 

sectional normal force coefficient 

C N 

total normal force coefficient 

Cn 

yawing moment coefficient 

C P 

pressure coefficient 

c t 

tip chord 

C T 

thrust coefficient 

c y 

side force coefficient 


*>Co 

= — r- , per radian 

Cn 6 

- , per radian 

Cyg 

~iCy 

= — -=r- , per radian 

C £p 

* ■»<*/»<&) 

^np 

- " >c " /»(■&) 

c yp 

» ■»<*/•(&) 

C&r 

- ^/*<&) 

C n r 

- 

Sr 

-■»<*/*<&) 

C t,tip 

tip suction coefficient 



f 


n 


G(x) 

l 

M 

N 

n 


N F 



P 

r 

R(x) 

s t 

tn* 


T 

u,v,w 

V 


strength of fuselage flow singularities of n th 
Fourier mode 

tip suction singularity parameter 
fuselage length 

Mach number, or number of integration stations 
number of integration points 

number of Fourier modes of induced velocities matched 
to satisfy the fuselage boundary condition (excluding 
the zero-order mode) . 

number of fuselage control or integration stations. 

the induced normal velocity of the n c h Fourier mode at 
the fuselage control station i due to a unit strength of 
flow singularity at the integration station k. 

the induced normal velocity at the control point l due to 
a unit strength of the horseshoe vortex density at k. 

roll rate, rad. /sec. 

radial coordinate, or yaw rate in rad. /sec. 
fuselage radius 
tip suction per unit length 
reference wing area 

the induced tangential velocity of the nth Fourier mode 
at the control point i due to a unit strength of flow 
singularity at the integration station k. 

■ 

nondimensional induced velocity components in the x, y 
and z directions, respectively. 

velocity 
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w 


x,y,z 


*P 



z cCx) 

a 

0 

6 

Y 

X 

A 

r 

r(x,y) 

X 

f 

0 


wing section width 

rectangular coordinate system with positive x-axis 
pointing downstream, positive y-axis pointing to the 
right and positive z-axis pointing upward 

the fuselage station behind which the potential flow 

ceases to exist and is taken to be 0.378 + 0.527 *■£/£ from nose. 

the fuselage station at which the rate of change of the 
fuselage cross-section area first reaches its maximum 
negative value. 

caunber coordinate 

angle of attack in degrees 

= V< - Moa 


sideslip angle, radian 
vortex density 
taper ratio 
sweep angle 
dihedral angle 
accumulated circulation 
= V~/Vj 
= yU COSC* 

nondimens ional velocity potential 
density 

angle between camber line and chord, including flap 
angle, if any. 


4 



Subscripts 


a 

c 

f 

ff 

fj 

fw 

£ 

3 

jj 

JJ 

JW 

Jf 

°j 

r 

s 

t 

w 

r 

8 


oo 


antisymmetric 

chordwise 

fuselage 

fuselage control stations being influenced by fuselage 
flow singularities 

fuselage control stations being influenced by jet 
vortices 

fuselage control stations being influenced by wing 
vortices 

leading edge 

jet flow 

jet flow perturbation due to jet vortices 

jet control points being influenced by jet vortices 

jet control points being influenced by wing vortices 

jet control points being influenced by fuselage flow 
singularities . 

external flow perturbation due to jet vortices 
radial direction in the cylindrical coordinate system 
spanwise, or symmetric 
trailing edge, or tail 
wing 

circulation 

circumf erential direction in the cylindrical coordinate 
system 

freestream 



2. Introduction 


It has been shown that upper-surface-blowing (USB) swept-wing 
configurations have unsatisfactory Dutch roll characteristics with- 
out stability augmentation (Ref.l). In an effort to evaluate the 
jet-interaction effects on the lateral-directional stability derivatives, 
and hence the lateral-directional dynamic characteristics, the present 
investigation was initiated. 

It is well known that some lateral-directional stability derivatives, 
in particular the dihedral effect (C^g) , depend significantly on the 
wing-body interaction (Ref. 2, Chapter 4). In addition, the fuselage 
will contribute to the derivatives Cyg, C n g, Cy r , and C nr . Therefore, 
the computer program developed earlier without the fuselage effect 
(Ref. 3) must be revised. 

The existing methods for representing the fuselage effect are 
based on either the image method (See Ref. 4) , or the surface singularity 
distribution (Refs. 4,5,6). The former method is applicable only to 
long cylindrical bodies with constant radius. On the other hand, the 
latter method can deal with bodies of arbitrary shape. However, the 
computer time and the memory requirement would be greatly increased. In 
order to keep the computer memory requirement to a minimum and yet have 
the capability of dealing with arbitrary bodies of revolution, a method 
based on the axial distribution of flow singularities will be developed 
in Section 3.1. The fuselage effect on the wing-jet interaction will 
be formulated in Section 3.2. Many ideas used in computing the lateral- 
directional stability derivatives in subsonic flow have been discussed 
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in the NACA and NASA literature, mostly based on the lifting-line 
theory. Some of these ideas will be reformulated for the use in the 
present lifting-surface theory (Ref. 7) in Section 3.3. 

3. Theoretical Formulation 

The basic assumption in the formulation is that the linear, 
potential, subsonic flow theory is applicable. The quasi-vortex- 
lattice method of Reference 7 is used to formulate the wing and jet 
boundary conditions (Refs. 8 and 9). 

3.1 Fuselage effect without jet interaction 

It is assumed that the fuselage is a body of revolution with 
arbitrary radius distribution. It is not restricted to a slender body. 
The fuselage camber is not accounted for in the present program, although 
this restriction can be easily removed. 

To represent the fuselage effect, the axial distribution of G.N. 
Ward's vortex multiplets (Ref. 10) will be used. According to Reference 
10, the velocity potential due to the axial distribution of vortex 
multiplets with strength f*(^) is given by 


r)«- 


I 

47r 





T 


(l) 


where cos n 0 is to be used for the symmetrical cases and sin n 0 
for the antisymmetrical cases, such as in computing the lateral- 
directional stability derivatives. The coordinate system for Eq. (1) 
is indicated in Figure 1. Note that if n=o in Eq. (1), the resulting 
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z 



Figure 1. Fuselage Cylindrical Coordinate System 


velocity potential is for the familiar line source distribution. Before 
Eq. (1) is differentiated to obtain the velocity components, it is 
advantageous to integrate the integral by parts to obtain 



The boundary conditions in symmetrical flow can be written as 


follows: 


On the wing surface: 


'Iz ' "*z 


sino( 


( 4 ) 
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On the fuselage surface: (see Ref. 11, Chapter 9) 


Z^£ + = coso< - sin* c«s0 (5) 

”■> r *r a* 

Note that both <f> v and <f>£ in Eqs. (4) and (5) are nondimensionalized 
with respect to V w . In Eq. (5) , is needed. It is obtained from 
Eq. (2) by differentiation to give 

= ~4?r ¥ ^ sn$ \ ^ .«*) 

J (6) 

Eq. (6) can be simplified as follows. By the following transformation. 





1 - C03 


»,) 


(7) 


the integral becomes 




where 



ill 

r r> V 


n & o 


( 8 ) 

(9) 


“ ^f (x ' r ^}t) (10) 

25± = 2ik r 

*» r "ir V *; J* ) (U) 


/ 
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( 12 ) 


and the integral was reduced to a finite sum through the midpoint 
trapezoidal rule„f n (xf t ) in Eq. (6) can be expressed in terms of 
f,i (O as follows. 




4 .1 


a I * 


It follows that Eq. (6) becomes 


H = _j_ 

*»r -+?r 




+ *f*~*i* Jr 




The fuselage surface boundary condition, Eq. (5), will be 
satisfied at Np control stations specified in a similar way as in 
selecting the wing control points (Ref. 7): 


^ ~ e,. = ^L i<:= i,...,k/ F 

3 4> 

► '4" J i f'/iei nta ac! T'V. n t 


Note that — 11 is expressed in a Fourier cosine series. The wing 
'ay- 

induced normal velocity, , will also be Fourier-analyzed. Let 

^ Y 

be a matrix with elements being the induced normal velocity 
of the n fc h Fourier mode at control station c due to a unit f' n k* 
Similar definition can also be applied to the matrix CA/j-yJ] due to 
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the wing effect. The detailed expressions are described in Appendix 
A. After Eq. (14) and the Fourier-analyzed are substituted into 
the fuselage boundary condition Eq. (5), and equating each Fourier 
mode, the following system of equations can be obtained: 


tX’Hf.'! -I (16> 


etc. 

In the above equations, all equations have been multiplied, through 
by R n * , where R(x) is the fuselage radius. This is numerically 
expedient in view of the fact that^fy^r , which appears in C » 

contains r*’ +l in the denominator. G n in the above equations comes 
from the first term in Eq. (14), and is defined as 


** ~ ^*$*^*3 $.****) (19 

Note that in each of the equations (16)-(18), there are Np equations 

t 

because of Np control stations, but involving (Np+ 1) unknowns, 
with k=l, ..., Np and f n (xf£) . However, f n (xf£) can be obtained 
directly by evaluating Eqs. (16)-(18) at x=x^£, i.e., the fuselage 
nose. The results can be shown to be (see Appendix B) 



( 20 ) 


■^V =o- 




1+3, 


_■ ,, ->A \ o> 

*xt + z: »W 


t 


)' ' 

**\'*''\ 


( 21 ) 



s. 



o 


7) Z.Z 


(22) 


where 


f, = 


1 I (23) 

V If 'X s *#* 

Note that in applications, it is not necessary to compute g^ because 


if 75 AE. J _o , it implies -^v'C^) * s finite so that g ( is nonzero. 
Ix-Xte s : ... 

■“ f i A 

However, ^ n (x^=o whatever g ( is. On the other hand, if is 

finite, this implies « if R(*jg) = o, and hence = 0. Eq. 

(16) can be solved independently of other equations. 

The wing surface boundary condition, Eq. (4), can be written in 
terms of the influence coefficient matrices, r*.wj. r< 3 . tviVj. 
as follows : 

t»U? \ M iI^-3 H'i + 

“ * • " r 1 • •-** • > •• “• ■ 

= \ 1— - “ H, 7 C 3 {$o\\ m] 

’ * 

where is the contribution of and its detail expression can be 

found in Appendix A. Eqs. (17), (18), and (24) can be written in the 
format of an augmented matrix for solution. Assuming n=2, it follows 
that 


4K 


etc. 


(^*w3. 3 £*443 

U) 

frCl LTC) ( °3 

1 *' 1 

k^3Lo3 Ck 1 ^] 

X' 


Ih 

'*X 


-smo<- ] 




(Sym ) 1 


J 


(25) 
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Note that in Eq. (25), G, contains / w through Eqs. (19) and (21). 
That particular term should be combined with the appropriate term 
on the left hand side before the solution is attempted. In all 
configurations tested so far, the contribution of to the final 
solution of / w has been found to be relatively small. 

For the antisymmetrical case due to sideslipping, steady 
rolling and steady yawing, the boundary condition can be formulated 
in a similar manner. Consider the wing boundary condition first. 
Referring to Fig. (2), the induced normal velocity, which must be 
produced to cancel the upwash due to dihedral in a sideslip, is 
given by 


_ 

( 26 ) 

where it is assumed that the wing boundary condition is to be 
satisfied on the mean surface only and T, the dihedral angle, is 
small enough so that the normal unit vector to the wing surface is 
essentially in the z-direction. In steady yawing, the sidewash 
produced on the wing surface would vary 






Figure 2 Normal Velocity Components on Wing-Body Surface due 
to Lateral-Directional modes of motion 
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from location to location and is given by rx as shown in Fig. 2. 
Therefore, it can be considered as in a "variable sideslip" situation 
with the following boundary condition to be satisfied: 

« = ' ( ' ' t) sln r = (ffeX^r) s ' m r <27 > 

Similarly, the upwash on the right wing in a positive steady roll is 
given by py. Hence, the wing boundary condition is 

2i =-±3- ,2SS 

-»Z V. '•2-V.) k <28) 

On the fuselage surface, the sideslip boundary condition is given 
by (see Fig. 2) 

~ — (~~p> Sin 6) ( 29 ) 

Similarly, in steady yawing, the expression would be 

(30) 

On the other hand, in steady rolling, no normal velocity will be 
generated by p-motion. Hence, 

2±. = 0 (31) 

In the case of antisymmetrical flow, the cos n 6 in Eq. (14) is to be 
replaced by sin n 0 and the wing-induced velocity on the fuselage 
surface will be developed in a Fourier sine series. The augmented 
matrix in this case can be written as, again for n=2. 


rtj„j z ki^] 


u, i 

'•■fw. 

,(?•>- 


X> - 


tu'gw to] 




-£ 5i *r+(^(^) sl < 




(32) 


^ o 

where the corresponding G n for this case are all assumed zero for 
simplicity. For S’ -derivatives, both p and r are assumed zero in 
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Eq. (32) and j3" is set to 5°. For p- or r- derivatives, either Ik 


r L 

or — — • is set to 1 and other flow parameters are set to zero. 




3.2 Fuselage effect with Jet Interaction 

In the wing-jet interaction theory described in Refs. 8 and 
12, the solution is obtained by solving the jet-off case first and 
then the additional effect due to jet interaction. The same concept 
applies here. The jet-off case has been described in Section 3.1. 
The jet effect of the symmetrical case is obtained by solving the 
following system of equations (Ref. 12) for n=2 : 

Pressure continuity on jet surface 
£ ” fC'"') C “T( < 4') > '[ $3 w3(o)$ 

- { f.'i = \ T 

Flow tangency on jet surface 

Flow tangency on wing surface 


Flow tangency on fuselage- first Fourier mode matching 

u = {-**«, s 

Flow tangency on fuselage- second Fourier mode matching 


(35) 


(36) 


+lT ' ,l « y ^\ - 1 »i (37) 

The velocity potential in Eqs. (33) and (34) is for the wing-body 
combinations without the jet effect. D]^ in Eq. (33) and P^ in Eq.(34) 
are the contribution from f^(xf£) and the detail expressions can be found 
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in Appendix A. 

In antisymmetrical case, Eq. (33) remains unchanged in form, 
except that all influence coefficient matrices must account for anti- 


symmetry 


and 4 > is interpreted as the antisymmetrical velocity 


potential for the wing-body combination without the jet effect. 


Furthermore, the right hand side of Eqs. (35) and (36) will be re- 


placed by zero and that of Eq. (34) by the following expression: 

RHS of Eq. (34) for antisymmetrical case =-^-/3 ccsSj + 

+ nr) sin (38) 

where 0j is defined in Fig. 3 and the effect of Mach number non- 



Figure 3 Jet Surface Boundary Condition in Lateral 
Mode of Motion 


3.3 Forces and Moments 

3.3.1 Symmetrical aerodynamic characteristics 

Once the vortex density distribution on the wing is obtained, the 
wing aerodynamic characteristics can be calculated by integration. The 
procedure has been reported earlier (Ref. 8) and therefore, will not be 
repeated here. The only difference is that the computed vortex density 
is referred to V* here, instead of Vo, cos a in Ref. 8» 

The pressure distribution on a body of revolution may be computed 
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by the following formula (Ref. 11). 


C Kf) = ' 2U '"[< | - M -) u '’'+ V v +w*-] (39) 

where u, v, and w are nondimens ional perturbed velocity components 

referred to the freestream velocity. In the fuselage cylindrical 
coordinate system, 


= % + v r v 

(40) 

where, referring to Fig. 1, 


! '”>4* 

= y ~ - s‘ino4 sin © 

(41) 

UJ. = ^-pr + Smo( cos© =. ««*<< - 3 -^- 

(42) 


with the first term coming from the combined effect of the fuselage 
and the wing. The last relation in Eq. (42) was obtained from the 
boundary condition, Eq. (5). Therefore, Eq. (39) can be written as 


1 ~ C 1-HOM- ( si* ©J" (43) 

The velocity component u due to the fuselage itself is formulated in 
Appendix C. Once the pressure distribution C p (f) is obtained, the 
fuselage sectional normal force coefficient referring to the radius 
can be calculated by integration along the circumference (see Fig. 4): 

= - 1 !, sw-* (44> 

One alternative in making the approximate integration in Eq. (44) is 
to divide (o,tt) into two arcs and integrate the results seperately: 


Eq. (45) is useful if the wing is at any arbitrary location on the 
fuselage. The wing-fuselage junction will be the dividing point. 



Figure 4 Decomposition of Pressure Force on the Fuselage 



With the sectional normal force coefficient calculated, the total 
normal force coefficient is given by 

where Xp is the fuselage station behind which the potential flow does 


( 46 ) 


not exist, and is related to the fuselage station x'p at which the 
rate of change of the fuselage cross-section area first reaches its 
maximum negative value (Ref. 13, also List of Symbols). To integrate 
Eq. (46) numerically, it is first transformed to an angular integration: 


* = + 3* ( i - eos<f>) 


4 


/* _ C * rc , sih 


(47) 


(48) 


2 

With (xf£-Xf t )rcN^^ (sin<j>)/2 expressed in Fourier cosine series, Eq. 
(48) can be exactly integrated. This is developed in detail in Appendix 
C. Similarly, the moment coefficient can be obtained from 


«L 


^Cf)* “H* 


(49) 


The relation between Xp and the fuselage geometry has been discussed in 
Ref. 13 and an empirical formula has been given. The same method is 
used in the U. S. Air Force Datcom. It should be noted that Eqs. (48) 
and (49) do not include the second order effect in the angle of attack, 
or any vortex lift effect. 

As shown in Ref. 11, the first harmonics in the last term of Eq. 
(43) will not contribute to the lift. Since the higher harmonics in 
jL should be small in magnitude, they are ignored in computing the 
lift and moment in the present computer program. 
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3.3.2 Lateral-Directional Stability Derivatives 

The theoretical prediction of lateral-directional stability 
derivatives has been discussed previously by Queijo (Ref. 14) with 
Weissinger's lifting-line formulation. More recently, Hancock and 
Garner have discussed in detail the mathematical theory and applications 
for computing the second-order forces and moments with the lifting- 
surface theory of the kernel function type (Ref. 15). The vortex-lattice 
method (Refs. 5,16, 17) and the pressure panel method (Ref. 6) can also 
be used for this purpose. In the present method, the quasi-vortex- 
lattice method (QVLM) of Ref. 7 will be used to predict the edge forces. 
It is assumed that the spanwise vortex density tfy has been predicted 
in accordance with the QVLM. 

(a) Streamwise vortex density distribution / x 

The basis of the present method is the calculation of the 
streamwise vortex density ( }f x ) distribution on the wing. According 
to the conservation of vorticity. 


2JL . 

'tx 

By integration. 


'z/u 




(50) 


Jim j£c*«,)22l 

•w * ^ 


(51) 


****' (52) 
If Eq. (50) is regarded as a first order differential equation for 
/x> then the second term on the right hand side of Eq. (51) is the 
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initial condition obtained from the fact that at the leading edge, 
the vorticity vector should be parallel to the leading edge. Let 





It follows that 


(53) 


V * 

* “* 'i'l 

3 (54) 

The problem of determining / x consists of two steps: Firstly, using 

the predicted /y-valu es at the vortex locations, calculate F(x,y) 
values at the same locations. Secondly, differentiate T(x,y) numerically 
to find / x . It should be remembered that jfy-values are defined at only 
a finite number of locations in the chordwise direction. To use these 
finite number of values in the integral Eq. (53), it is best to develop 
/y sin 0 in a cosine Fourier series, where 0 is the angular coordinate 
in the chordwise coordinate transformation used in the QVLM: 

x =■ x z + (i - cos©) (55) 

Let 

r / *1 (56) 

-He) = x, sine = + 5: ^ 

The Fourier coefficients can be computed as 


- ij « 7 > 

a i = ± , j* i <58) 



H-i 

IN 


nr 


(59) 
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Using Eqs. (55) and (56), Eq. (53) becomes 


r ^, V =-¥ ) 5 = ( 60 ) 

Eq. (60) has been found to converge very well with respect to N, the 
number of Yy~values defined in the chordwise direction. 

With discrete values of r(x; c ,yj c ) defined over the planform, 
theoretically it is possible to pass an interpolation surface through 
these "points" , with (x^y^, T^) defining a point in the three-dimensional 
space. At each point, a tangent plane to the interpolation surface 
exists. This tangent plane can be constructed if the slopes in two 
directions (to be specified) at a point are known. These two directions 
can be conveniently chosen to be x- and y'-axes (see Fig. 5), where 



x 

Figure 5 Coordinate Transformation 

the y'-axis is along the straight line on which the predicted Yy-values 
are situated. It follows that 


X' = X Cos /\ — ^si*A 

(61a) 

#* = X 5i«A + J tosA 

(61b) 

X = x'ccsA + a' si* A 

(62a) 

3 = “ x'sinA-f Cos A 

(62b) 
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where A is the sweep angle of the y'-axis. Therefore 


•?>r __ 2H 2. 2JEL _ 


'ix 


I 




- S '*A~ + ^A^, 




( 63 ) 


dp _ -*r 22! 4- 2 H , * a 

^,X TX' ^>X ^ y 'Zx n>x' + 1H 

T> P 

Substituting from Eq. (64) into Eq. (63) results in 

~t>X 

2L. = _-u A 2T _L_2H 

^ T>x + CO sf\ -*y 


(64) 


(65) 


-ft p ^ P 

Since — — = -y v (x,y), or y v can be obtained through Eq. (65) if 

-ox ” '"»H x 

can be computed. In the present method, has been obtained by a 
trigonometric interpolation formula developed in Appendix D. 

(b) Tip Suction 

The tip suction can be predicted by using the trigonometric 
interpolation formula for developed above. The tip suction per 

w 

unit length of the tip chord is given by (Ref. 18) 

s *00 = TrSSoo (66) 

where G(x) is the singularity parameter of T(y) at the tip and is 
related to f(y) by 


GCx) J!*! /'--ifc 2 


(67) 

The detail of computing G(x) by using the trigonometric interpolation 
formula developed in Appendix D is formulated in Appendix E. 

(c) p -derivatives 

The distribution of lifting pressure coefficients can be written 
in terms of the wing upper surface perturbed velocity as (Ref. 19) 




t 


"»X' 


(68) 
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where x' is measured along the freestream in sideslip (see Fig. 5). 

For small sideslip angle g, Eq. (68) can be written as, by using Eqs. 

(62) with A replaced by g. 

The change in the lifting pressure in sideslip results from the wing- 
body interaction with or without dihedral (see the boundary conditions 
(26) and (29)) in sideslip and the term -gy x in Eq. (69), where y x is 
from the symmetrical loading. The second term can also be explained 
by Kutta-Joukowski theorem in that the sideslipping velocity will 
interact with the streamwise vortices (y x ) to produce the positive 
lifting pressure on the right wing in positive lift. On the left wing, 
the lifting pressure is negative, thus creating a rolling moment. Hence, 
the rolling moment coefficient due to sideslip ("dihedral effect") can be 
calculated as 


ri^hr wirv^ 

where y x is regarded positive if as a vector it is printing downstream, 
■^Yy/dp is due to the wing-body interaction with or without dihedral and 
0 the camber angle (including flaps) for proper force resolution. 


To compute the wing contribution to Cyg and C n g, the change in the 
wing edge forces with the sideslip must be found. For a wing with di- 
hedral, or a wing-body combination with or without dihedral, an anti- 
symmetrical pressure loading will be induced. This will change the 
wing edge forces. The sectional leading-edge thrust coefficient in 
combined symmetrical and antisymmetrical loadings can be written as 



(Ref. 7): 


**, = -*■ 


2c/r-nr — rr (<4 + 0' 


(71) 


where C s is the leading-edge singularity parameter in the symmetrical 
loading and the subscript "a" denotes the antisymmetrical case. It 
follows that 

Tr / r 

S ~ " “v'i- nu (72) 

where the subscript 8 denotes differentiation with respect to the 
sideslip angle. C a can be computed in exactly the same way as in computing 
C s . With Ac X g obtained /iCyg is given by 


- 4C *a Vn Ajj 


(73) 


The total contribtuion to Cyg and C n g can be computed by integration 
In addition, the incremental pressure loading due to sideslip will produce 
antisymmetrical change in the tip suction. According to Eq. (66), the 
tip suction can be written as 

S* * 'frS ( 5s ± 

where the plus sign is for the right wing; the subscript s denotes the 
symmetrical loading and a for the antisymmetrical loading. It follows 
that 


*fcs 


=^5S s (±^) 


(74) 




where — - can be obtained in the same manner as G_ . Therefore, the 

-ip ^ 

wing side force due to sideslip will be contributed from the following 
sources : 

(1) Contribution from the change in the leading-edge suction 


** v« 3M' wmo 

This is shown in Fig. 6(a) 


( 75 ) 
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(2) Contribution from the incremental pressure force 

This is illustrated in Fig. 6(b), where it shows that 
with the lifting pressure acting normal to the planform, 
there will be a component of the pressure force contributing 
to the negative side force. 

(3) Contribution from the induced drag (page 14-3, Ref. 20) 

The induced drag is assumed to act in the direction of 
the freestream with sideslip. This is illustrated in Fig. 

6(c). Hence, if Cj>^ is the induced drag due to the symmetrical 
loading, then 

= - C Vl /£ 

(76) 

Of course, the viscous drag will also contribute to Cyg. But 
this is neglected here. 

(4) Contribution from the Tip Suction 


This is given by Eq. (74) 



Figure 6 Wing Contribution to Side Force 
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Similarly, the wing contribution to the yawing moment due to sideslip 
results from the following sources: 

(1) Contribution from the change in the leading-edge suction. 

This can be seen from Fig. 6(a). By taking moment of the 
forces Ac x g and Acyg about the z-axis, AC n g can be obtained. 

(2) Contribution from the incremental pressure force 

The incremental lifting pressure force, being assumed to 
act normal to the camber surface, will be resolved in a direction 
parallel to the chord plane. After multiplying by the moment 
arm y and integrating, the contribution to C n g can be obtained. 

In addition, the side force due to the dihedral shown in 
Fig. 6(b) will also contribute to C n g. 

(3) Contribution from the tip suction given by Eq. (74) 

The fuselage contribution to the side force and yawing moment can be 
calculated in the same manner as in computing Cjj(f) and C m (£) due to 
fuselage. In the present case, the sectional side force coefficient is 
given by (see Fig. 4). 


-*7r 


c y(-f)"“ r i < “hf/ ra0)sme 

O 

Therefore, Cy(g) and C n (f) are computed as 

^m) = S rCy w 


u 




W X 


Since the fuselage surface pressure Cp^ is computed at 6=5°, 


(77) 

(78) 
(79) 

the results 
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of Eqs. (78) and (79) will be divided by 0 to obtain the fuselage 
contribution to Cyg and C n g. 

(d) p-derivatives 

The roll damping derivative Cj^p is computed by integrating 
the antisyinmetrical pressure force induced by the roll rate multiplied by 
the spanwise moment arm, assuming that py/2\/ m =\ • The pressure force 
distribution is antisymmetrical, because with the positive roll rate, the 
right wing will see incremental upwash, while the left wing will be subject 
to incremental downwash. This can also be seen from the boundary condition 
Eq. (28), where y is positive for the right wing and negative for the 
left wing. 

The side force and yawing moment due to roll rate, Cyp and C n p, 
respectively, are assumed to be due to the pressure distribution on the 
wing surface only. In other words, the fuselage contribution is ignored. 

Of course, the wing-body interaction has been included. The contribution 
to Cyp and C n p results from the following sources: 

(1) Contribution from the change in the leading-edge suction 

The incremental sectional leading-edge thrust coefficient 


is given by Eq. (72) with appropriate changes: 


AC *r ** - t A- K 




-s 

Cos Ab 


where Ac xp and C a are the derivatives with respect to ( Ib/zVc*) • 


Again, ACy* is obtained from 


^C /(J = -KnAj (81) 

With A c x p and A Cyp known, their contribution to Cyp and C n p can be 
obtained by following the same procedure described previously for 
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Cyg and C n g . 


(2) Contribution from the tip suction 

The idea expressed in Eq. (74) is still applicable 
for this case. 

(3) Contribution from the incremental pressure distribution 

This contribution is non-zero if the dihedral angle is 
not zero (to both C_ and C„ ) and if the camber is 

yp n p 

present (to C n p). The concept described in relation to Cyg 
and C n g is still applicable. 

(e) r-derivatives 

The incremental pressure distribution due to yaw rate consists 
of three components. 

(1) Due to yawing, a backwash ry is produced (see Fig. 2). This 
will interact with the symmetrical y y to reduce the lifting pressure 
on the right wing: 


*<r = - m / 
r r v. s 


t> h 


rk/iv, = I 

(2) The sidewash rx will interact with the symmetrical y x to 


(82) 


produce the following ACp r : 


AC = - i 
V V y 


- -12LV (83) 

L> 

! rb 4v 00 =i 

(3) Wing-body interaction with or without dihedral. See the 


boundary conditions (27) and (30). 

With the incremental antisymmetrical pressure distribution obtained, the 
wing rolling moment due to yawing can be computed immediately. 

The calculation of C yr and C nr follows the same procedure of 
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computing Cyg and C n g by recognizing the fact that a wing in yawing 
will produce the "variable sideslip" effect, because the sidewash 
varies on the wing and the fuselage. 


4. Numerical Results and Discussions 
4 . 1 Fuselage alone 

In computing the fuselage surface pressure, Eq. (43) is used 
with the following modification. Since the theory is a linear one, 
near the nose the pressure would be unrealistic. This can be seen from 

J t? 

the second term in Eq. (43), where could be infinitely large near 

the nose. Therefore, this term is modified to be 

4R 


4R "<4 X 

cos o( -777 Replaced by (cos<x + u) 


(84) 


** s 'VhT§? 

This modification is similar to the Riegels factor used in the airfoil 
theory (Ref. 21). Figs 7-9 show the camparison of the pressure distribution. 
Except very near the nose, the present method appears to be reasonable 
in predicting the pressure distribution as compared with the experiment 
(Ref. 22). It should be noted that without the modification indicated 
in Eq. (84), the pressure within about 10% of the nose would be more 
negative than that presented in the figures. 

4.2 Wing or Wing-body characteristics without jet interaction 

It is well-known that the fuselage interaction will increase the 
loading on the wing. Fig. 10 shows the comparison of the predicted 
incremental loading with experimental results given in Reference 23. 

The present results were computed by including up to second-order 
Fourier modes in satisfying the fuselage boundary condition (i.e., 
n=2 in Eq. 6). The agreement is good except near the wing-body juncture. 
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Experiment. . 

Exact Potential Flow 
. [...Theory .7 

Linear Slender Body 
Theory 

Present 


Figure 7. Comparison of Experimental and Theoretical 
Pressures on a Prolate Spheroid of Fineness 
Ratio 6 at a * 0 and M =* 0. 


Figure 8. Comparison of Experimental and Theoretical 
Pressures on a Prolate Spheroid of Fineness 
Ratio of 6 at a = 5.6° and M_ =0. 0 = 0°. 
















Before comparing the prediction of lateral-directional stability 
derivatives, it is necessary to make certain that the tip suction 
prediction is accurate, as the Cy - and C n - derivatives depend on it. 

For this purpose, a cropped delta wing of aspect ratio 1.333 and taper 
ratio 0.5 is chosen. The present method is formulated in Appendix E. 

The present results are compared with Reference 18 in Table 1. It is 
seen that the agreement is good and the present method yields reasonably 
stable value with respect to the vortex arrangement. 

With the accuracy of the tip suction prediction established, it is 
now possible to compare some of the predicted lateral-directional stability 

Table 1. Comparison of Predicted Tip Suction Coefficient for 
a cropped Delta Wing of AR = 1.333 and A =0.5. 


a = 1 radian 

N c x N a Lamar Present 

1.400 

5 x 20 1.361 

7 x 20 1.377 

9 x 20 1.387 

5 x 12 1.364 

7 x 12 1.384 

9 x 12 1.395 


derivatives. In Table 2, the p-derivatives for three wings are presented. 
The present results are seen to agree quite well with Gamer's (Ref. 15). 
In Table 3, the predicted derivatives, C^ T and C^g, are compared with 
Queijo's method. Again, the agreement is reasonably good. 

Finally, the results for a swept wing— body combination of Ref. 25 
are presented in Fig. 11. xi^which is related to Xp in Eq. (46), is 
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taken to be 1. It is seen that C £ g is well predicted, except at 
higher angles of attack. On the other hand, for C ng » only the trend, 
not the absolute level, of the high- 


Table 2. Comparison of Predicted p-Derivatives with Gamer's 
Calculation (Ref. 15). Based on Body Axes. 


Wing 



Garner 

Present 



v° 

c £p C np/CL 

Cyp/a C £ p C 

np/ c L 

Rectangular, 

AR=4,Moo=0 

1.374 

-0.3360 

-0.168 

1.391 -0.3367 

-0.171 

Rectangular, 
AR=4 , Mo>=0 .866 

1.945 

-0.3794 

-0.140 

1.963 -0.3802 

-0.145 

Tapered-Swept , 






AR=2,A £ =60° 
Moo=0 . 7806 

2.733 

-0.1854 

-0.365 

2.6027 -0.1848 

-0.336 

Table 3. 

Comparison of Predicted Derivatives 

with Queijo's 



method. 

Based on 

Body Axes. 



Derivatives 


C 

)ueijo 

Present 


V/CL 


Ref. 14 
0.4 

Ref. 24 

0.384 


C£B/ C L 



-0.47 

-0.485 



wing configuration is predicted. For the low-wing configuration, the 
predicted results at low angles of attack are more negative than the 

experimental values. The C vo derivative, not shown, follow the same 

” p 

trend of C^g. In other words, for the high wing configuration, the 
trend of Cyg is well predicted, but the predicted level is only about 
half of the experimental values. The discrepancy is most likely due to 
the effect of nonlinear lift and moment on the fuselage as discussed in 
Reference 13. The difference in C y g and C R g for high and low wing 
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configurations may have resulted from the wake displacement and de- 
formation. However, this can not be proved in the present investigation. 

4.3 Effect of Jet Interaction 

The experimental results on the lateral-directional stability 
derivatives with jet effect are scarce. When available, they are mostly 
for the complete configuration (Refs. 26, 27). In Ref. 26, some results 
with T-tail off are presented. However, the fuselage has its after body 
upswept. To compare with this set of data, it is assumed in the present 
theory that the fuselage is a body of revolution, with elliptical fore- 
body, cylindrical section within the wing root chord and parabolic after 
body. The maximum diameter is obtained by measurement of the graph, as 
the dimension of the fuselage is not given. From the model static test, 
the thrust efficiency can be determined to be 87%, and the jet deflection 
angle to be 44° at <5f=60°. to predict C L with the jet off as closely 
to the experimental values as possible at low angles of attack, the 
actual flap angle used in the program was assumed to be the same as the 
jet deflection angle, i.e., 44°. The predicted Cl is obtained by adding 
the computed incremental Cl due to jet effects to the experimental jet- 
off values. The results for the lift curves are shown in Fig. 12. It 
is seen that the lift is slightly overpredicted at C^=3.6 and under- 
predicted at C<£=1.8 for the assumed configuration. The sideslip 
derivatives for the same configuration are presented in Fig. 13. The 
main difficulty in predicting Cyg and C n g in the present case is that 
because of the upsweep of the fuselage after body, the flow separation 
would occur earlier than the case with a body of revolution. In the 
present computation, x^jlis taken to be 0.8. The agreement is seen 
not as good as that shown on Fig. 11, in particular, for the dihedral 
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Experiment 

O c T =o 

□ C t = 1.8 
O C x =3.6 

Present Theory 



o 

o 


I I l 

6 8 10 
g* 

Figure 12. Comparison of Predicted Lift Curves 

with Experiment for an USB Configuration 









effect. The exact reason for the disagreement is not known. From the 
theoretical point of view, for this high-wing configuration can be 
positive only if the symmetrical lift is negative. However, some ex- 
perimental values are positive. In addition, the experimental results 
with the jet on are not consistent in that at a=0°, with the power 
reduced, the value of does not come closer to the jet-off values. 

It seems that some nonlinear viscous effect may have been present under 
the high lift conditions. The theory has a tendency to overpredict 

at high lift, as can be seen from Fig. 11. It should be noted that 
the direct power effect due to the momentum transfer at the engine 
inlets on Cyg and C n ^ is not included in the figure. When it is included, 
both Cyg and C n ^ will be more negative. However, the uncertain contribution 
from the nacelles and fuselage will make this refinement unnecessary. 

As indicated in Ref. 1, the swept wing configurations with the USB 
jet have unsatisfactory Dutch roll characteristics. Therefore, it is of 
interest to see how the leading-edge sweep will affect the lateral- 
directional stability derivatives. For this purpose, the unswept high- 
wing configuration of Ref. 26 is again used, with the leading-edge sweep 
angle increased to 30° without changing the aspect ratio. The results are 
plotted against the circulation lift coefficient in Figs. 14-16. Only 
those derivatives which show significant changes are presented. A general 
observation of the results is that although the circulation lift coefficient 
is greatly increased with the jet on, those lift-dependent derivatives are 
not proportionally changed for this configuration. It is seen that the 
swept wing configuration has better directional stability (more positive 
C n g) with or without the jet on than the unswept one, and has large 
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dihedral effect. The sweep largely increases Cyp, mainly because of 
the contribution from the leading-edge suction. C n p is more negative 
at a given Cu 1 . From Fig. 16 y it is seen that the leading-edge sweep 
provides better yaw damping (more negative C nr ) , and has large C£ r . 

Large positve Cjj, r is detrimental to the spiral mode, but more negative 
C nr is favorable. From the above results, it is seen that the 
unsatisfactory Dutch roll mode of the swept wing configuration is 
caused by the large dihedral effect without power. With the power on, 
the dihedral effect is further increased and makes the situation worse. 

5. Concluding Remarks. 

A jet interaction theory for calculating the lateral-directional 
stability derivatives for wing-body combinations has been developed. 

The fuselage effect is represented by the axial distribution of G.N. 

Ward's vortex multiplets. The predicted results show good agreement 
with experiment and other available theoretical methods for configurations 
without the jet effect. Because of lack of data for comparison, the 
accuracy of the theory with the blowing jet effect cannot be established. 
Further comparison with available data in the future would be needed if 
the refinement of the theory is to be made. 
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Appendix A 


Influence Coefficient Matrices for Wing-Body-Jet Interaction 

(1.) Effect of Fuselage on Wing 

Consider a wing position on the fuselage as shown in Fig. 

A. 1. At any point on the wing, the upwash induced by the fuselage 
can be shown to be 

VTs: - V & si*© (A.l) 



Wing 


Figure A.l Induced Velocity on Wing due to Fuselage 


Using the formulation with G.N. Ward's vortex multiplets, T/ r is 
given by Eq. (14) „ _ _ 

r -jr ■** L sIhn©J l 

+ H f (A . 2) 

where X^.- > and c<?sn© is to be used for symmetrical cases 

and si**© for antisymmetric al cases. To find V a , Eq.(2) is differentiated 

with respect to © . 

i -ai , , f -nsinhei r 

v «= =-fv r ,c OTAS ii tF * w ' r '*«)-R.Kn* H )3f h («: f .p 

(A - 3> 
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■*F« 


In the above equations, and are defined as follows 


F„u,r, 5 )=- 1 . 

v*n y 


(A. 4) 


r rV^-pN^r' 

Substituting Eqs. (A. 2) and (A. 3) into Eq. (A.l), the upwash 
expression due to each Fourier component can be obtained. For 
example, in the case of symmetrical loading, the upwash due to 
Fourier component n is 


(A. 5) 


* r - s? { f c*v 




where it is understood that =z —■ (x, r, » etc - 


(A. 6) 


The influence coefficient matrix £ is then defined as 


x *>r 

w f ^ 


w&J «•*»*) JL > -l 

** 1 **,1.9*3 a A/p T>r VJ 

+ s>n6 i V">t»»nne A -? i) -v . r 

cr.^- r.c v3 


(A. 7) 
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The term ^j ( appearing in Eq. (24) is given by the remaining terms 
in Eq. (A. 6) with : 



(A. 8) 


Figure A. 2 Induced Velocity on Jet due to Fuselage 

Referring to Fig. A. 2 where *n is the unit outward normal vector to 
the jet surface, it can be seen that 

ri s» j cos 8* + -ft »•■»* 8* (A. 9) 

The velocity vector induced by the presence of the fuselage can be 
expressed in terms of VJ. and as 



It follows that the normal velocity component induced on the jet is 
given by 


■u;= vti =(v r si*e + v e «>s$)«>s9'+(v r cose-v 6 sme)s'me , (A11) 


With Eqs.. (A. 2) and (A. 3) substituted into Eq. (A. 11), the influence 
coefficient matrix ^A/^T^can be obtained: 

l*“l ■+»Uinhe c /* s, "^Uf ( V-^(S|)] 

The term ?, in Eq. (34) is given by 

- «**.>- rscv3^ v -%• 

+ ^iniv „ce, n $J (A.13) 

To find the u- influence coefficient matrix £ (2) 

differentiated with X to give 

-+r 

= -±z( 

sl*i>9 





where 


If* 

~bX 




-L r , >0 

■/(*-'$?+ ft* 


(A. 15) 


"a x **«) ^ K > r > Xft ) 


(A. 16) 


It follows that the influence coefficient matrix T^j-p J given by 

3h»>3 (*.«,■ 


and in Eq. (33) is 



(2.) Effect of Wing-Jet on Fuselage 


(A^ 18) 


The general expression for the induced velocity at any point 
due to the wing-jet vortices has been given before (Ref. 8). Consider 
a unit horseshoe vortex on the wing or the jet surface. At a 

fixed number of points to be specified on the fuselage circumference, 
both the upwash wand the sidewash V will be computed. The normal 
velocity at a given point at station t will be given by (see Fig. A. 3) 


V* =• W-cosS + V5in0 


(A. 19) 



Wing 


Figure A. 3 Induced Velocity on Fuselage due to Wing-Jet Vortices 
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Let 1^, be developed in a cosine or sine Fourier series without the 
zero order term, depending on whether it is a symmetrical or anti- 
symmetrical case: 


v; 


» 21 a J 

»=i L 


c cshO 

a n ] 

Then, the Fourier coefficients will be given by 


(A. 20) 


2 if'* C ) 

= . Ue 

n o i s»»>,e j 


(A. 21) 


To illustrate the numerical integration, consider the integral 
= ^ l^cosn&clS 

C 

The interval ( O ,*)r ) will be divided into two parts: (0,6) and 
(8,,7r) (see Fig. A. 3). Hence, 

2 . c r' 1r 

a »= j J "u; oosn&<A0 + } cos *8 

^ 0 J 


" * l \ f “ sn8 « + ^ f V nj } 


(A. 22) 


where the integrals have been reduced to finite sums through the 
midpoint trapezoidal rule and 


e ‘ = ^ 6 ' 


> & — | 


(A. 23) 


9 = S,+ l ^OCTT-g i) ,j =1 


N/j 




(A. 24) 
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The angular positions indicated by and are the locations at 
which the induced velocities are to be computed. With a n obtained, 
the influence coefficient matrices I and are then given 

by 



(A. 25) 


(3.) Effect on Fuselage due to Fuselage 

The influence coefficient matrix can be obtained from 

Eq. (A. 2) as follows: 




(A. 26) 


where ^ {.x^) ’ t * le fuselage radius at *» . 



Appendix B 


Derivation of in the Formulation for the Fuselage Effect 

(1.) The Expression for 

The zero-order Fourier component of the boundary condition 
Eq. (5) can be written as 


^4 dft 

1 s tOSO< — — 

1) r 

From Eq. (6), it is obtained that 


(B.l) 




= «** 


(B.2) 


where the equation has been multiplied through with r=R(£. Now, 

Eq. (B.2) will be evaluated at According to Eq. (6) for the 

■SF 

expression —la. , the following limiting values can be derived, assum- 
ing R(*#)=0 » 


r ?hf*> =, 


r £ c V = 










r-±\ix) 




(B.3) 


(B.4) 


r~W = 7 > 

,r yc*-5? +/ iv- *-*»*, 


(B.5) 
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It follows that 


X=> X 




~4Trf f \ \ 

I * fl 

^ '"t &l) == COSot 

Therefore, 

,, in | 

r f s_ +'* c * s * * ** !*=*« 

**> i+T — 

(2) -f, (.***) 

The first-order Fourier component of Eq. (5) is 


-a4> o; 'a4> <0 

+ = - sino< 

ir 'if 


(B.6) 


(B. 7) 


(B. 8) 


Again, Eq. (6) gives 


“ts 


~J»W 

+ 

* '»r 


= -R smof, 


r*Roo 

(B.9) 


r=ROe> 


To find , Eq. (B.9) will be differentiated once with X before 

the limiting values as *— » a:^ are evaluated. Hence, 






r0c-A^f-f^r'] 3/x i 


V 

r-*R(x) 

*-* x H 


l-l + [-l + >} 


= 0 


(B. 10) 
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where is assumed small compared with ( ) as 

&>$<*>> • + 7f=^c + ( *'V? 55^ 

- w^+A-^l j 


»l +S+ y 9 

' '" ri»«jsf^“ 1,1 

■tt 


(B.ll) 


Similar to Eq. (B.10), it can be shown that 


gZr'fcv"} 


(B. 12) 


It follows that 

I 


•+7r \ ^ x fi) \ ^ ^ K> 


x=x. 


tt 




« 


- _1R TGTL_ 


si>io< 


x=x. 


■fc 


*.i v = - r si ”« + * C/«» 1 ) 


(B.13) 


By inspection, it can be determined that with jet on, Eq. (B.13) will 
be modified to be 



(B. 14) 
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Similar procedure can also be followed to derive the expressions 
for > 71 5- 2. • However, the algebra will be more complicated. 

On the other hand, by observing Eq. (B.13)^tcan be determined that 
),nj z , would depend on the higher Fourier components of 
the wing or jet induced velocities at the pointed nose. Clearly, 
these higher Fourier components are zero. Hence, 


■£>(*#) = ° , *2 2 . (B . 15) 

Eq. (B.15) can be easily shown by the slender body theory. 

It should be noted that will affect mainly the pressure 

distribution near the nose of the fuselage. Because its effect is 
small on the over-all characteristics, it will be assumed zero in 
the prediction of lateral-directional stability derivatives. 
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Appendix C 


Calculation of Fuselage Aerodynamic Characteristics 

(1.) Calculation of the u-perturbed velocity on the fuselage surface 
due to the vortex multiplets. 

As will be shown later in this Appendix, the sectional normal 
force coefficient is to be computed at the integration stations 
(i.e. , at the stations where -f^J-values have been computed). 
Therefore, dj, and u are also to be computed at the same stations. 
For better accuracy, u = - ■ * ■ will be computed as follows. 

Differentiating Eq. (2) for the fuselage velocity potential with 
X , it is obtained that 



The integral will be integrated in the following way. Consider 7)= 0 . 
Then 


^C5) = - 1 


(C. 2) 


f - - — 1 . t j 1'cp-f.w Jt , rtt 


(C.3) 
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It follows that 


/ 

£ V f -°?- 

T $ VCx-^f+^V 1 




(C. 4) 


Note that x is the location where _J5is to be computed and in the 


present case coincides with ^ 


Similarly, for different n , it can be shown that 


£ , / tosft6\r 

~ 4?r ( Sb, * e ) y * ^ wT 

^ _ ^bs) 4 ?} 

& **« 


(C.5) 


^ l ^ I 1 (c * 6) 


2Ev^ „ _ J_ j>-18 r 

** r" yr«-«***‘r* 


r /C*~5) t +/sV‘ 

r^t 

3 


(C.7) 
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(2.) Fuselage normal force and moment coefficients 

Once the pressure distribution on the fuselage surface is 
computed, it can be integrated to give the sectional normal 
force coefficient . Integrating C M along the fuselage 


axis will produce the total normal force coefficient. 

Accord- 

ing to 

Eq. (48), it is given by 



%>= VS rc *« EmHf 

(C. 9) 

Let 

fc40 = 

(C.10) 

^^»)can be developed in a cosine Fourier series: 



-R4»)= a Q + £ a. co s j<f> 

(C.ll) 

where 

o F ■*=! 

(C.12) 



(C.13) 


A 2-& - I - . 

” «/ f ^ 

(C. 14) 

Once the Fourier coefficients are obtained, direct integration 

of the 

series will give : 



, A Wp . 

= J fCfMf sin 

(C.15) 


Similarly, the moment coefficient can be obtained as 
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Appendix D 


A Trigonometric Interpolation Formula for the Calculation of 
Streamwise Vortex Density 

As shown in Eq. (65) , the streamwise vortex density can 

-ji- 
be computed if — is known on the wing surface, where is 

a 

the inclined distance as shown in Fig. D.l. Assume that 



the r -values defined in the section of width W are to be used 

■y 

will be used: 


for interpolation to find . The following two transformations 



(D.l) 


*'= J£(I-W S 0)+ a' (D.2) 

Eq. (D.l) transforms the interval into ‘^(o^Tr) and Eq. 

(D.2) transforms the inclined interval ^■W*) into 



Since it is well known that p possesses square-root singularities 
at both tips, it is numerically more accurate to develop p sin <£> , 
instead of P itself, into a cosine Fourier series. 


= F s '*4> 


(D.3) 


+ X- a i tosjO 

J=l 

where the Fourier coefficients are given by 


(D.4) 




(D.5) 


a r** 


o ^ ' 


(D.6) 


where the conventional, not midpoint, trapezoidal rule has been 
used because the p -values are computed at the locations which 
are not appropriate for the midpoint rule. is defined as 


e * = 


(D.7) 


If Eqs. (D.5) and (D.6) are substituted into Eq. (D.4), it is 


obtained that 


M M-l 


Kfc) = XI casj ft, c«S 3 '0 

1 3=1 -*=l * 

+ + e M ) 

Differentiation of with respect to $ gives 

2 M M-l 

~ M s,H j © 


(D.8) 


(D.9) 
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Eq. (D.9) is to be evaluated at 0^ , where 


A 


©. = ±2L 

M 


Using known trigonometric relations, Eq. (D.9) becomes 


4^. I M-l M 

7 *. mT?, 2*£*£*W< VM] 

~ m Z. a 

5"“ I 3"* 3 


However, it is known that (Ref. 21) 


M 

ZL i 

3=1 


M s'm(M-H ) x - (M+l)sinM* 

2.(*os x - |) 


Hence , 


M 

£ * E *'*i C G>$ + SO - S'nj (f>* ~ ©*■ >] 



[ 


4&i 



^(efc+eo-i 

cos (- ^^ -) 
si-n( 


sm(8t-^) 
tosC^- ©_.)-! 



= 'MC-0 




sih 9^- 


, -«r+-C 


cos©,, - cos 0^. 
Similarly, it can be shown that 


J! < . , _ M sinCM-vO©i ■" 

2- 1 S '»,J 9i = — 

1=1 z ( cos ©^ - l) 


M(-l)*~ s'mSt’ 

z C l - cos 


(D.10) 


-«<)) 

(D.ll) 


(D.12) 


(D. 13) 


(D. 14) 



rftx 


3 [ s ' m 3 (&*,+ ©;) “ s '»i 


_ Ml-0 y sin (6 M t6^) _ »m( 6h-fc) ~i 

1 L Ms (e M +e v )-i «»(e M -Si)-l J 

= Mt-0 M+i ±*i- 
\ + ccs e - 

If k=i, Eq. (D. 13) will be revised to be 


M 


2: isins-je^ = 

3~l a(cos^-i) 


M sinxS^ 
a(l- c«s2©-) 


M_ u>s 0 t ' 

2- si* 0 * 

Therefore, Eq. (D.ll) becomes 




M-l 


d $. t*(-0 


•Vk si*0: 


c-os0<; 


■& = I 


+ — 5. — — 

c**^-**^ * i 3i *e* 


+ ±5 c ~ oSiw ^ _j., , 

I-CCS0. I + cos 0. 

where the prime on the summation sign implies that the 
is to be omitted. 

Differentiating Eq. (D.3) with ^ gives 

_ -»P . L . - . 

• 7 ; — Sin <P t r 00 -S <P 

T ‘ 


p 

Solving for — ■*- , it is obtained that 


(D.15) 


(D. 16) 


(D- 17) 
term with 


(D.18) 
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2JL = J_ 5 il_ ££< 

^ 3 'i. <13 


_ . 4^ 

? - r 608 1; -JY 009 ‘ 


= _l_Sil _i X«£ii-« s Al 

s'm^l 4^ WVt^B^ -L sj*^ J 


(D.19) 


Eqs. (D.17) and (D.19) represent the desired formula for computing 2J1. 


It should be noted that ^ and ^ can be obtained by interpolation. 
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Appendix E 


Calculation of Tip Suction 

According to Eq. (66) , the tip suction per unit length of the 
tip chord is given by 


S*<*) = Jrs 6 X C*> (E.i) 

where g(x) is related to the circulation through the following 
relation (Ref. 18) 


Near the right tip, Eq. (E.2) implies that 

ru, j) = - 46 w/E /,_ ( 4 .) 

1 7 

Eq. (D.19) shows that 


. i *»P _ 4 3 a A 

sm«j> — — • = — — 1- cos A 

4© IV'smB -k. 

But at the tip. 




_ -ax 2r 2E_ a 

> -a-Y T* H cos/ * 


►V •»* ’’S 




y 

because ;r— 'sr-jt, vanishes at the tip. Hence, the objective is 
to find the limiting value of the following equation: 


(E.2) 


(E. 3) 


(E.4) 


(E.5) 


5 . , 2. Pmsj ) 

~ *f> wsme 3 


(E.6) 


where W ** W'm>*A has been used. Differentiating Eq. (E.3) with 
respect to ^ gives 

_ 2.G I 


■* s ^ /^¥) 


(E.7) 



It follows that 


lim sin<J> 




(E.8) 


U 4^1,1, -^V^H-0 
-»i **•.+ W'-c-fey yr 


(E.9) 


where 


si>i ^s/ ^^ ^see Eq. (D.l)). Furthermore, Eq. (D.17) 


gives 


il 1 

<i ©«, W si* 


n s<H 


^ Cos 8-- 0)30- ^ 


+ J_ p s ' m< k 

1 u/ 1 « r» 


l-tos^- 


L r ♦ r C-0 




(E.10) 


where the second term on the right hand side of Eq. (D.17) has 
been combined into the first term because i will not be equal to 
in the present case. Note that at the tip i = M and Q i =r qy 
(see Eq. (D.10)). In addition. 


**** = »- 


(E.ll) 


1+ cos ft = li>M 3 -(\ — 


(E.12) 


where Eq. (E.ll) is obtained by multiplying Eq. (D.2) by cos A • 
The last term in Eq. (E.10) has the following limiting value as 


I*,* r < _ \\ r * ~ 4g v / T/*' 

y-*i. l+cosp- 2. ( I— 

3- 2. ( I - jtlls ) 

— = 1 ^ 


(E. 13) 
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where the third step was obtained by applying L’Hopital’s rule. 
Eq. (E.10) now becomes 


*3 


V*+M 


— — = rr i n sin 4* V 4ip s :a H) M 

Wsj n ^ w -*=} * H tos^ + | W *• *t„ 


'• 2 


+ 


/F 


(E.14) 


Substituting Eqs. (E.8), (E.9) and (E.14) into Eq. (E.6) gives 


yz 

a 


^£'nsm£ c '° 


-few 


* * l+cos©^ W e 2. 


■ ft - 1 


From Eq. (E.15), G(x) can be determined. 


The local tip suction coefficient is defined as 


(E.15) 


_ 3*<X> _ 2/V5 x 


%-kip j_ 


xS\Lst 


> v.«i 


(E .16) 


The total tip suction force is found from 


r *»* V 

F * “ 2 J '1rSG*-Ax sa c x } A © 




Sc ,%! ¥ p.‘ s H 

■fc=i 


It follows that the tip suction coefficient is 


(E.17) 


F* 


r = — -z- 
+i r -L?v 2 




2-C* M ^ . 
* «■ ? s « 


* =l (E.18) 

Eq. (E.18) can be used to estimate the side edge contribution to 

the vortex lift as shown in Eef. 18. 



